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Abstract 

Fixed energy inverse scattering theory has been used to define central and 
spin-orbit Schrodinger potentials for the scattering of 5 eV polarized electrons 
from Xe atoms. The results are typical for a range of such data; including 
energies above threshold when the potentials become complex. The phase 
shifts obtained from an analysis of the measured differential cross section and 
analyzing power has been used as input data. Both semi-classical (WKB) 
and fully quantal inversion methods have been used to extract central and 
spin-orbit interactions. The analysis shows that information additional to the 
set of input phase shifts extracted from this (and similar) data may be needed 
to ascertain physical potentials. 
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I. INTRODUCTION 



A knowledge of the interaction between colliding quantum systems is central in many 
applications of scattering and has relevance for use in other diverse fields of study. Such inter- 
actions have been sought in various ways with the method of numerical inversion common. 
In the numerical inversion approach, the parameter values of a purely phenomenological 
parametric form, chosen a priori to be the (central, local) interaction between the colliding 
entities, are determined by variation until a best fit to measured data is found. Global inverse 
scattering methods fl|] form an alternative procedural class with which to analyze the same 
data. With global inverse scattering methods, the interaction between the colliding pair 
is extracted from the data without a priori assumptions about the shape of the potential, 
although it may belong to a certain broad class, and the validity of the dynamical equation 
of motion (the Schrodinger equation) is assumed. Potentials so obtained we define here- 
after as inversion potentials. Application of various global inverse scattering methods have 
been made in the past for electron-atom |2| , for atom-atom || , and for electron-molecule Q 
systems, but none of those methods permitted extraction of spin-orbit effects. Recent de- 
velopments [§-0] have provided means by which spin effects can be treated. In this paper 
we present and describe results for central and spin-orbit potentials obtained by global in- 
verse scattering methods. In particular we consider an approach based upon the WKB 
approximation method |j]|§ and one based upon the Newton-Sabatier (N-S) theory Jl|. 

The data of interest come from the very high quality crossed beams experiment of Gibson 
et al. ||. In that study ||, a phase shift analysis was also made and the phase shifts so 
obtained were the input quantities to our studies. Those phase shifts are purely real in 
line with the unitarity constraint with the energy below the first threshold. Thus we have 
obtained purely real inversion potentials. Extension of the approach to deal with energies 
above threshold and, concomitantly with complex potentials, is straightforward. The key 
feature about inversion potentials, given numerical accuracy in calculations and stability of 
the solution, is that when used in Schrodinger equations they lead to the same phase shifts 
as are input. 

The data set we have chosen to use is intriguing for a number of reasons. First, the 
differences between the values of 8f and 5f (the ± superscripts denoting j — I ± ~) ex- 
tracted by the phase shift analysis || are not large indicating that the spin-orbit interaction 
is not strong. As a result the approximation method of Leeb, Huber and Fiedeldey || 
(LHF hereafter) can be used with confidence. Indeed the LHF scheme is accurate through 
second order in Born approximation and has worked well in some nuclear scattering data 
analyses [|K|] where the spin-orbit effect is much stronger than in the case we study. With 
the LHF approximation, both the N-S and the semi-classical WKB methods of inverse scat- 
tering theory can be used to specify the electron-Xenon (e-Xe) potentials. Exact quantal 
inversion methods to get the spin-orbit interaction are known J^J^], but with this data the 
LHF approximation should be adequate and the inversion process is facilitated by its use. 
Second, the phase shifts of significance are not many in number and so this may be another 
case with which phase shift values at unphysical rational values of angular momentum are 
required in the inversion process to achieve a stable result |TT| . A third reason for interest is 
that the s- and p-wave phase shifts from the analysis || of the scattering data have negative 
values. All other (physical) phase shifts of significance are positive quantities. As phase 
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shifts are ambiguous to modulo n, an equivalent completely positive valued and monoton- 
ically decreasing phase shift set can be formed by the addition of 2tt and 7r to the s- and 
p-wave values respectively. With either of the original or the modulated (integer I) phase 
shift sets as input, the N-S inversion method per se gives the same inversion potential. How- 
ever, that potential has a short ranged repulsion; such being required (12| to give significant 
negative s- and p-wave phase shifts. The WKB inversion method on the other hand, does 
discriminate between these sets since to specify the WKB inversion potential interpolated 
functions, <5;(A), are required as input. Of course if the N-S method is extended to use phase 
shifts at rational values of the angular momentum, found for example by interpolation of the 
original sets and of the those with s- and p-wave values adjusted by 2ti and tc respectively, 
then the inversion potentials from the two cases must differ. 

Thus we give two new elements in this paper. First we deduce by inverse scattering 
theories, central and spin-orbit potentials for 5 eV electron-Xenon atom scattering. Second, 
we show that information additional to the physical phase shifts (i.e. those determined by 
usual phase shift analyses of scattering data) is needed to identify the the most likely physical 
inversion potential. In the next section we give a precis of the LHF approximation for phase 
shifts as well as of the inverse scattering methods used, the N-S and the WKB semi-classical 
schemes specifically. Then in Sec. Ill, we discuss the origins and characteristics of the 5 eV 
e-Xe scattering phase shifts that have been used as input to our inversion studies. The e-Xe 
potentials that result are presented and discussed in Sec. IV and we draw conclusions in 
Sec. V. 



II. FIXED ENERGY INVERSION METHODS 



In this section we give brief outlines of the methods used in the calculations, the results of 
which we report later. First we set out the LHF scheme by which the spin-orbit interaction 
can be defined from two independent (spin-less) inversion calculations. This identifies not 
only the special phase shift sets {Si} and {Si} of the method, but also defines the central 
and spin-orbit potentials in terms of the results of inversions of those new phase shift sets, 
V and V respectively. Then we give the salient features of the Newton-Sabatier and the 
semi-classical WKB inverse scattering theories which we have used to specify the (V and V) 
potentials. 



A. The LHF approximation 



While exact quantal inverse scattering theories that yield central and spin-orbit inter- 
actions from input scattering phase shift sets exist @,0], Leeb, Huber, and Fiedeldey || 
developed an approximation scheme to transform the input phase shift sets so that more 
facile quantal inverse scattering methods, such as the N-S scheme and the semi-classical 
WKB approximation ||, can be used to give results from which central and spin-orbit po- 
tentials can be extracted. Those more facile schemes do not allow for an angular momentum 
dependence in the intrinsic equation of motion, such as given by a spin-orbit potential. They 
are designed only to provide central local potential functions. 
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The LHF method is based on the assumption that the contribution of the spin-orbit 
potential to the phase shifts can be evaluated using a Distorted Wave Born Approximation. 
This technique has been formulated specifically for spin ~ particles incident on spin zero 
targets and is accurate to second order in the Born expansion ||. The approximation 
identifies first the special phase shift sets {Si} and {Si}, and then defines the central and 
spin-orbit potentials in terms of the results, V and V respectively, from inversion of those 
new phase shift sets. 

For spin | particles incident on a spin target, and allowing central and spin-orbit 
Schrodinger potentials, the scattering is defined by reduced radial Schrodinger equations 



d 2 1(1 + 1) 1 
dp 2 p 2 E 

where, with p = kr, 



V em (p)+atV,M i>f(p)=o, (1) 
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If the spin-orbit term is relatively weak, the usual scattering phase shifts, Sf, can be 
expanded in powers of af. Specifically || 

6? = Sr + afC^ik) + {4fC?\k) + .... (3) 

While the leading term, Sf en , is due solely to the central component of the potential, V cen , 
higher terms must be considered to define the spin-orbit properties. But the specific analytic 
forms of the coefficients, C l - , do not have to be known to extract the central and spin-orbit 
potential values. The LHF approximation is initiated by considering combinations of S + and 
S~ from which separate inversion potentials can be estimated. The relevant combinations 
are 

~ Sl = oTTT { (/ + 1)<5 < + + U i) = 6r + /(/ + l)C ? ] + • • • ' (4) 



and 



k = — {ist + + i)s r } = sr - + (i 2 + i + \)c\ 2) + ... (5) 

To first order in af, these new phase shifts and their inversion potentials, V and V, are 

Si = Sr ~ V ~ V cen (6) 
Si = Sr - Cf } ^ V ~ V cen - ~V S0 . (7) 

As these new sets of phase shifts can be inverted independently using any of the conventional 
techniques, the central and spin-orbit components can be identified then by 

V cen (r) « V(r) (8) 
V so (r)*2[V(r)-V(r)} . (9) 
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B. The N-S method 



Since the Newton-Sabatier inverse scattering theory and applications have been widely 
reported, only pertinent points of the scheme are presented herein. A full treatment on the 
development of this method is given elsewhere 

The N-S method is one of the most successful of the fixed energy inversion methods. Very 
recently, it has been applied successfully to electron-helium atom scattering Jl3[] using as 
input, experimental phase shifts of Nesbet [|T?[ at low £- values and dipole polarization phase 



shifts at the higher ^-values. But it is known [15| that fixed energy inverse scattering theory 



requires the S matrix (equivalently the phase shifts) as a function of the angular momentum 
variable if one is to define uniquely the scattering potential. This equates to knowing the 
S matrix exactly at all of the (infinite) set of physical /-values as then the unit step in the 
quantum number is infinitesimal against the range. Most studies of the fixed energy inverse 
scattering problem, and notably those involving the N-S method, have been applied using 
only the values of the S matrix specified at a finite set {l m ax) of physical angular momentum 
values. In cases where there are relatively few important partial wave phase shift values to 
be used, it may be necessary to extend [TT^] the usual N-S formulation to include rational 
values of angular momenta to form the matrix inherent in the N-S scheme [|l]. 

However it serves to consider for this section, just the integer values of angular momentum 
for which the Schrodinger equations take the form (p = kr being dimensionless) 



d u ( P ) ^r(p) = Ki + 1) wo*) (io) 



where the operator 



D u (p) = P 2 



(11) 



with U(p) = V(p)/E cm , where E cm = (hk) / (2/j) and p = kr. The solutions are subject to 
boundary conditions 

flip)— t A sinfp-^7r + 7 A (12) 



with 7; being the relevant phase shifts to be taken as input quantities. The N-S method 
gives as output 

U(p) = U Q {p)--~K(p,p) (13) 
pdpp 

wherein Uq is a reference potential and K (p, p) is the Jost transformation kernel which can 
be written as the infinite sum of solution function products, 

K(p, P f ) = J2c l fl(p)fl°(p) (14) 
l 

The solution functions (to D u ) can be expressed by the Newton equations [|TJ 

0r(p) = «p)-E c <'Mp) flip) as) 
i' 
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where 



MP)= [ P <P?°(p'W(p)^dp> (16) 
jo p 

These equations are of central importance. From them one can determine the unknown 
quantities, Ai and (qAi), by matching asymptotically to the defined boundary condition 
solutions for p > p , p being the value at which the unknown quantal interactions are 
presumed to be vanishingly small. There is also a presumption that the solution functions 
of the reference potential are completely known so that the initiating L-matrices can be 
defined exactly. The reference solutions are obtained from 

D"°(p)tf°(p) = l(l + l)ff°(p) (17) 

with 

0r°(p) p - co sin(p-i/vr + ^) (18) 

where 5f is a reference input phase shift. With the normalization and expansion coefficients 
so given, the complete solution functions can be determined from Eq. (|T5|) at all p < p$. 
Thereby one gets the Jost transformation kernels and thence the sought after potential. 



C. The WKB method 



In the WKB approximation with A = I + ~, scattering phase shifts are defined [|lT,|l"2| by 

71 f°° 

5(X) = -X-kr + \K x {r') - k\ dr' (19) 

L Jr 



where K\(r) describes the local momentum through the interaction region and ro is the 
classical turning point. Thus the scattering potential is embedded in K\(r') and inversion 
amounts to an integral transformation. To effect such a transformation it is convenient to 
consider the deflection function 



9(A) 



d5{X) 
! dX 



(20) 



where now A is taken as the angular momentum variable. This deflection function satisfies 
an Abel-like equation, found by applying the Sabatier transformation, 



a = kr 



1 - 



Vwkb{t) 
E 



(21) 



to Eq. (pD. One finds 



5(X) = -±r -M oda 
K J 2EJx Ja* - X 2 



lx V° 2 - A 2 
where Q(a) is a quasi-potential defined by 



(22) 
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Q{a) = 2Eln {^) • (23) 

The Abel-like integral equation for S(X) can be inverted to give 

Q{?) = — t / rr=m XdX > 24 

7r a da Jcr Vcr 2 — A 2 
which can be written in terms of the deflection function as 

Provided there is a one to one mapping of the transcendental equation 

fcr = aexp (^r)' (26) 

and the energy i? exceeds that at which 'orbiting' occurs, i.e. 

E > V(r) + \r^- , (27) 
2 dr 

then the Sabatier transformation equation provides the relationship from which the scatter- 
ing potential can be found, namely 



V W kb{t) = E 1 1 - exp 



(28) 



For large a, the quasi-potentials decrease so that with a — > kr, Q(a) — > V{r). As 
a — > however, the quasi-potentials diverge and the transforms then lead to the lower limits 
r — > tq (the turning point radius), V(r) — > E . However, in practical cases the validity of 
the WKB approximation breaks down at a radius larger than r , when the transcendental 
relationship between a and r becomes ambiguous. 



III. SPECIFICATION OF SETS OF PHASE SHIFTS 

The 5 eV e-Xe phase shifts determined by Gibson et.al. have interesting structure 
notably that while the s- and p-wave phase shifts are negative, for all other Z-values they 
are positive. The filled circles in Fig. |l| depict the the phase shifts that have been extracted 
from the data. The top most graph identifies the phase shifts associated with the j = I + \ 
angular momentum set while the bottom panel contains those associated with j = I — |. 
It is evident from the data displayed in this figure that the there is only a small difference 
between the 5f sets. The largest difference occurs with the p-wave phase shifts, and that is 
only of order 0.1 c . 

As the phase shift analyses of the e-Xe scattering data gave negative values for the s- 
and p-wave phase shifts, one can expect scattering potentials that have a short ranged 



repulsion. But for e-Xe scattering it is known that the potential should be attractive at all 
radii and especially so near the origin where the incoming electron should feel essentially 
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only the presence of the nucleus. Then one would expect the phase shifts for low /-values 
to be positive. Such can be formed with the phase shift values having a monotonic decrease 
with I by the addition of 2ir to 5q and it to 5f. We define such modulated values as the 
7r-adjusted phase shifts hereafter. Naturally multiples of any integer amount may be used, 
but this new set is the simplest. The new (7r-adjusted) values are shown by the open circles 
in Fig. |I| Associated with such phase shifts are purely attractive interactions which are 
expected in the physical potentials for electron-atom scattering. 

To investigate the effect of additional input on the form of the inversion potentials, the 
two data sets were interpolated. Several interpolations were made seeking suitable input for 
the different inversion schemes. A many point interpolation was made on each phase shift 
set to obtain the input for the WKB inversion scheme. Values of the phase shift functions 
had to be found at quite small step sizes, Al, since in the WKB we have to evaluate not 
only the deflection functions but also integrate over them (numerically). A step size Al of 
0.01 was used. Also two extended sets of input phase shift values were generated for use 
with the N-S scheme. One had Al = 0.5 and the other Al = 0.2. This was done to assess 
the effect of differing numbers of non-physical input on the N-S inversion potentials. The 
sets of interpolated phase shifts obtained using Al = 0.2 are displayed in Fig. [I], with the 
solid and dashed curves giving the results of those (spline) interpolations. Clearly the phase 
shift functions so specified are no longer equivalent and so we expect any inversion process 
that requires such functions as input to give different inversion potentials. 

IV. RESULTS AND DISCUSSION 

The results we have obtained using the N-S inverse scattering theory are discussed first 
and, subsequently, those found from our WKB study of the chosen two phase shift functions 
are considered. We present in three subsections, the potentials that result, the phase shifts 
found from solutions of the Schrodinger equations containing those potentials, and the cross 
sections that ensue in each case. 

A. The results from N-S inverse scattering theory 

Our N-S studies have lead to six inversion potentials; three found by using the original 
phase shift values of Gibson et al. || and the other three obtained by using the 7r-adjusted 
phase shifts. For each case, we first calculated the N-S inversion potentials using as input 
solely the phase shifts corresponding to the physical I angular momentum set I 6 0, 1, 2, 7. 
Such results we identify as case 1 results, e.g. case 1 inversion potentials from the 7r-adjusted 
phase shift sets. Two other calculations have been made with N-S inversion. First the N-S 
inverse scattering theory equations have been solved using the discretization Al = 0.5; the 
results we identify by the designation case 2. The third set of N-S calculations were made 
using Al = 0.2 to give what we term case 3 results. 
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1. The N-S inversion potentials 



Physical arguments dictate that the e-Xe scattering potential for 5 eV electrons should be 
real (the energy is below the first threshold) and to be attractive, with long range behavior 
of — |a|/r 4 and a short ranged one of —Ze 2 /r. One may also expect that the intermediate 
range potential would be essentially a monotonic function between the extremes as the 
charge density of the atom is believed to be a smooth function 

The potentials resulting from the inversion of the phase shifts based upon the original 
set || are all strongly repulsive at small radii and hence are not considered physically 
significant. They also have marked oscillation in both their central and spin-orbit results. 
When used in the Schrodinger equations however, the solutions do reproduce the input phase 
shifts quite well; comparable to the results we show subsequently. But as the inversion 
potentials are not consistent with the form of e-Xe potential dictated by knowledge of that 
scattering system, we consider those inversion results no further. 

In Fig. ^] the potentials obtained by inversion of the 7r-adjusted phase shift values are 
displayed. The top and bottom segments portray the central and spin-orbit components of 
the potentials respectively. The dashed, long dashed, and solid curves represent the case 
1, case 2, and case 3 potentials respectively. The dashed curves in this figure are identical 
to the results of inversion found using the original phase shift values of Gibson et al. 
That is as should be since within the N-S inversion scheme made with just the phase shift 
values specified at the physical angular momenta, modulo 7r adjustment means that one uses 
exactly the same expansion wave functions in defining the internal matrices. But the other 
cases have quite different outcomes. The potentials shown in Fig. § tend to the physical 
expectation and clearly demonstrate that the inclusion of greater numbers of phase shifts 
at non-integer angular momentum lead to smoother, more realistic potential forms. The 
inversion potential found using the phase shifts at solely the physical angular momenta 
does not represent a structure expected for 5 eV electrons on Xe atoms as it has the short 
ranged repulsion. However between 0.75 and 2.5 a.u., that (case 1 central) potential has an 
attractive well with a depth of -1.3 a.u. while beyond 2.5 a.u. it behaves approximately as 
r -4 . The spin-orbit component of the case 1 potential is very small, is weakly attractive in 
the vicinity of 1 a.u., mildly repulsive between 1 and 2.5 a.u., and after that it is essentially 
zero. 

Obviously the most realistic potential comes with the case 3 potentials found using the ir- 
adjusted phase shift sets. That concurs with the hypothesis that the inverse scattering theory 
result stabilizes with increase in the number of non-integer angular momenta phase shifts 
commensurate with numerical accuracy of evaluation. Essentially the angular momentum 
step size should be small in comparison to the number of significant partial wave input data 
i}max)- This case 3 potential, shown in Fig. || by the solid lines, has exactly the structure 
one would associate with e-Xe scattering. At small radii it is strongly attractive and of 
r _1 form, with a smooth transition to a long range r~ 4 tail. There is a smooth transition 
between those regions. The spin-orbit results also are more reasonable. The case 3 (with 
7r-adjusted phase shifts as input), spin-orbit potential is not as extensive as the others. But 
the spin-orbit potentials are all small in general (save for the naturally occurring divergence 
at the origin) and so these three results do not by themselves indicate convergence. The 
reproduction of phase shifts and observables, however, tend to suggest that the results we 
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show are reasonable. 



2. Reproduction of the phase shifts 

An indication of the success of an inversion scheme is to reproduce the input phase shifts 
from solutions of the Schrodinger equations using the inversion potentials. In this case 
such reproduction is reasonably good but not exact; perhaps being a measure of the LHF 
approximation. 

In Fig. H the original phase shift values |9[] are portrayed by the filled circles and are 
compared with those obtained (at integer angular momenta) using each of the three inversion 
potentials defined by the 7r-adjusted input sets. Case 1 results lie within the filled circles, 
case 2 values are shown by the open squares, and case 3 gave the results portrayed by the 
open circles. The three lines now are meant only to guide the eye by connecting the phase 
shift values arising from use of the three inversion potentials. 

The best reproduction of the phase shifts defined from experiment |§ is found by us- 
ing the case 1 inversion potential, notwithstanding that the potential contains unphysical 
characteristics. Essentially, the inversion scheme in this case produces a potential that has 
been defined from just the two sets of eight phase shift values found by Gibson et al. at the 
physical values of / = to / = 7. Case 2 and 3 potentials, on the other hand, were built using 
many more phase shifts specified at non-integer /-values and as the inversion potentials then 
seek to reproduce all of those extra values equally well, small variations in the results at the 
8 physical /-values from the sets of Gibson et al. can result. The choice of those additional 
(unphysical) phase shifts then is crucial if the resulting potentials are to reproduce the eight 
physical values very well. One has to balance the need for a sufficiently large basis so that 
the inversion potential has stabilized to the proper (physically credible) limit, against the 
numerical accuracy one needs to achieve with reproduction of the physical phase shifts and 
scattering data. 

A source of possible error in addition is the choice that must be taken for the phase 

shift at the unphysical point 5 J 2 . That value is needed in the calculations of both 5 and 
S and also in the N-S scheme. This choice has the potential to introduce error since its 
value influences the interpolation. A poor choice of this value can become evident when 
the inversion potential is used to recalculate the phase shifts, particularly for the low-/ 

o=- x j=+ 1 

partial waves. A reasonable choice seems to be to set the phase shift S 2 equal to 5 2 . 
Admittedly this is an arbitrary point. In this study, two choices were considered; the first 

being to take 5 2 = 5 2 and the other to use an Akima spline to determine the value by 
extrapolation. In this study allowing the Akima spline to determine this point was slightly 
more successful. 

3. The cross section from the N-S inversion potentials 

Although the potentials all look reasonably good, particularly that found from case 3 
with the 7r-adjusted input phase shifts, the further test of the inversion results is to see how 
accurately use of the potentials reproduce experimental data. This is displayed in Fig. £| 
wherein the experimental data (with error bars) as found by Gibson et al. [g] are compared 
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with the cross sections calculated from the three inversion potentials obtained using the ir- 
adjusted phase shift sets. Those results found using the case 1, case 2, and case 3 inversion 
potentials are depicted by the dashed, long dashed , and solid curves respectively. Further we 
display the results here on a linear scale to distinguish the bulk of the results in relationship 
to the error bars. Later when discussing the WKB calculated cross sections, we will display 
also the case 3 N-S comparison with data shown on a semi-logarithmic plot. That emphasizes 
the comparison of results with data at the large scattering angles and particularly in the 
vicinity of the minima as 120°. As one might expect from the reproduction of the phase 
shifts, we see in Fig. [| that a good reproduction of the cross-section data is found with the 
case 1 results. That result passes through the error bars of most data points. The case 

2 cross section has similar structure to the experimental data, but the shape is slightly at 
variance falling just outside the error bars associated with a number of the data. The case 

3 cross section however is in excellent agreement with most of the experimental data. In 
general it falls within most of the data error bars except at the larger scattering angles; the 
latter indicative of the phase shifts at integer values of I not being reproduced with sufficient 
accuracy. 

B. The results from semi-classical WKB inversion theory 

The WKB inversion results have been obtained by forming the S(X) and S(\) phase 
shift functions and using them to evaluate two quasi-potentials. The potentials that result 
are discussed in the first subsection. Subsequently we present and discuss the phase shift 
reproductions and the cross sections that result on using those inversion potentials. 

1. The WKB inversion potentials 

The inversion potentials we have found using the semi-classical WKB methods are pre- 
sented in Fig. |5|. Once again the central potential is shown in the top section of the figure; 
the spin-orbit potential in the bottom. Clearly two very different structures have been gen- 
erated. The result from using the 7r-adjusted functions has physically sensible characteristics 
but that found using the deflection function defined from the original phase shifts does not. 
Indeed the inversion procedure based upon the original phase shift set does not lead to a 
spin-orbit potential one can identify as anything sensible and so that is not displayed. 

The 7r-adjusted phase shift data sets give smooth monotonic phase shift functions and 
hence work well with the WKB scheme. Furthermore it leads to the form of potential 
(the continuous lines in Fig. ||]) one would expect for the central component of an e-Xe 
interaction, i.e. a smooth function which is highly attractive toward the origin and has an 
attractive r -4 long range behavior. The spin-orbit component of the potential would also be 
small and short ranged. Essentially the potential from this WKB analysis that is displayed 
in Fig. [| has this desired prescription. However, given that the input was not optimally 
smooth at large angular momenta (close inspection of the input data showed that there 
were small oscillations in the phase shift functions), there are small oscillations at large radii 
in the WKB inversion potential. Nevertheless, by smoothing, the long range WKB inversion 
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potential does behave on average like r -4 and the central component of this WKB inversion 
potential resembles that found using the N-S scheme. 

Only the central potential found from WKB inversion of the original phase shift data || 
is shown in Fig. [| (by the long dashed line). Quite evidently it is nonsensical. It exhibits 
ambiguous behavior at many radii, resulting from a loss of 1:1 correspondence between a 
and r in this case. Apparently if the input phase shift function contains a large degree of 
curvature then that input is unsuitable for use with the WKB procedure. This is emphasized 
when one considers the deflection function, quasi-potential, and the associated a vs r plots. 
The deflection function resulting from differentiating the phase shift function defined from 
interpolation of the original phase shift values, has a rapid variation as is evident from the 
top panel in Fig. ||. Consequently the quasi-potential will also be quite structured and that 
is shown in the middle panel of Fig. | As a result the stability condition breaks down at 
a fairly large radius. This condition, the correspondence between a and r, is displayed in 
the bottom section of Fig. [?]. Clearly between 3 and 5 a.u. there is an ambiguous relation 
with cr and so one finds ambiguous values of the associated WKB inversion potential in that 
region. Indeed one can only hope to ascribe a sensible result with this input to the WKB 
inversion for radii then in excess of 4 - 5 a.u. 

The realistic e-Xe potentials found by using the 7r-adjusted phase shift sets and with 
both the N-S (case 3) and WKB inversion schemes, are compared in Fig. [?[ The potentials 
found using the N-S scheme are shown by the solid curves while those found by WKB 
inversion are displayed by the filled circles. The latter are taken to the limit radius at which 
the 1:1 correspondence between a and r is maintained. Overall the (semi-classical) WKB 
potentials compare very well with those found using the (full quantal) N-S scheme. Clearly 
the central potentials look almost exactly the same. There are differences between those 
two which become apparent when the observables are calculated. On the other hand, the 
spin-orbit WKB potential is very smooth and weak in comparison to the corresponding N-S 
component. Also it is repulsive at all radii while the N-S potential has a small attractive 
well between approximately 0.3 and 0.6 a.u. 

2. Reproduction of the phase shifts 

Despite the very pleasing form generated for the potential using the WKB scheme, the 
reproduction of the phase shifts, although reasonable, is not as good as those generated 
using the N-S scheme. However this discrepancy is not disconcerting given that the WKB 
scheme is a semi-classical approximation. The phase shifts found from our WKB calculation 
are shown by the open circles in Fig. |8], and joining lines have been included solely to guide 
the eye. The data of Gibson et al. || again are represented by the filled black circles. As 
with the N-S case, part of the discrepancy between the calculated values and the data may 
be due to the ambiguity with the interpolation required to specify the phase shift functions 
and particularly for points below I = 1 for the j = I — | input set. 

A spline was used to determine the deflection function, and that influences the quasi- 
potential and ultimately also the potential. Given that points found from the j = I — \ 
input set are extrapolated for values of A < |, ambiguity in those values is inevitable. That 
ambiguity persists when these splined values are used to determine the input functions <5;(A) 
and 5i (A). That may result in a poor reproduction of the physical phase shifts, particularly 
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of the s- and p-wave values when the solutions of the Schrodinger equations specified with 
the relevant WKB potentials are found to complete the study loop. 

3. The cross section from the WKB inversion potentials 

Already, from the reproduction of the phase shifts, we suspect that with the WKB 
method, reproduction of the cross section will not be as good as that obtained using the 
relevant N-S inversion potentials. This is indeed true. The solid curve shown in Fig. |9] depicts 
the WKB cross section result which is compared therein against the Gibson et al. data || 
and also against the cross section determined by using the preferred case 3 N-S inversion 
potential (shown by the dashed curve). We now present the results in a semi-logarithmic 
plot since WKB cross section is similar in magnitude to the data, however it does not display 
quite the right structure at any scattering angle. The reproduction simply is not as good as 
that found using the N-S scheme. This figure also emphasizes the mismatch of the case 3 
N-S inversion results at large angles that we commented on earlier. 

V. CONCLUSIONS 

Inversion potentials for the 5 eV electron-Xenon atom interaction have been found using 
both the full quantal N-S and the semi-classical WKB inverse scattering theories. The 
results from application of the N-S inversion were very good. Several inputs were used in 
this approach, each containing a different number of phase shift values. From solutions 
of the Schrodinger equations specified with each of the potentials, 'inversion' phase shifts 
were extracted that reproduced well the starting (physical) phase shift values and also the 
empirical cross section. However, when the input was taken solely to be the (physical) phase 
shifts at integer-/ values of angular momenta, the inversion potential contains a marked 
repulsion at small radii. That is not physical. With increased numbers of input data 
specified at non-integer values of angular momentum, inversion produced a potential with 
sensible (physically credible) characteristics. 

As two disparate inversion methods find central e-Xe potentials that are essentially the 
same and have pertinent physical properties of the scattering system, we are confident that 
the N-S approach has (nearly) converged and that the central potential obtained by expanded 
7r-adjusted physical phase shift sets is the appropriate candidate for the local Schrodinger 
interaction. The spin-orbit potential is reasonable but more detailed investigations are 
needed before the characteristics found for it can be adopted with confidence. 

The similarity between the WKB and N-S inversion potentials also implies that the 
introduction of unphysical phase shifts in the N-S calculation is essential if a stable inversion 
potential is to be obtained in this case. That would be so with other energies measured || 
and, by implication, for any such electron-atom scattering at eV energies. There is the 
difficulty however, of accurately specifying the phase shifts for non-integer angular momenta. 
Obviously some kind of a priori information regarding the colliding system is necessary. In 
this instance simply 7r-adjusting the given phase shift data sufficed. 
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FIGURES 
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FIG. 1. The scattering phase shifts found from a phase shift analysis of 5 eV e-Xe scattering 
data. The filled circles are the results specified by Gibson et al. [Sj while the open circle are the 
7r-adjusted values. The results of interpolations of the basic two sets of (physical) phase shifts are 
portrayed by the solid and dashed curves. 
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FIG. 2. Potentials (central top; spin-orbit bottom) obtained from N-S inversion of the three 
sets of phase shift values described in the text and formed by using the 7r-adjusted phase shift 
values. The dashed, long dashed and solid curves depict the results designated cases 1, 2, and 3 in 
the text respectively. 
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FIG. 3. The 5 eV e-Xe phase shifts obtained from the potentials shown in Fig. || compared 
with the values specified by Gibson et al. fjj) . Note that the lines are simply to guide the eye and 
to identify the three cases. The notation is as used in Fig. ||. 
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FIG. 4. The 5 eV e-Xe differential cross sections obtained from the potentials shown in Fig. ^ 
compared with the data of Gibson et al. || (dots with error bars). The notation defining the 
results from the three cases is as in Fig. 0. 
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FIG. 5. Potentials (central top; spin-orbit bottom) obtained from WKB inversion of the phase 
shift functions described in the text. The solid curves depict the results found using the 7r-adjusted 
phase shifts to define the phase shift function; the long dashed curve gives the (central) potential 
found using actual Gibson et al. M tabled values for that purpose. 
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FIG. 6. The deflection function (top), the quasi-potential (middle), and the a vs r plot 
(bottom) from the WKB study framed upon the original phase shift values of Gibson et al. ||. 
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FIG. 7. Comparison of the WKB inversion potentials (filled circles) with those of the case 3 
N-S inversion study made using phase shift sets from interpolation of the 7r-adjusted phase shifts 
of Gibson et al. g. 
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FIG. 8. The 5 eV e-Xe phase shifts obtained from the (realistic) potential shown in Fig. || 
compared with the values specified by Gibson et al. Q. Note that the lines are simply to guide 
the eye. 
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FIG. 9. The 5 eV e-Xe differential cross section obtained from the WKB potential (solid 
curve) compared with the data of Gibson el al. |J] as well as with the case 3 N-S result (dashed 
curve). Both inversion studies used interpolations of the 7r-adjusted set of phase shifts as input 
data. 
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